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Multiple outflows in the planetary nebula NGC 6058 
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ABSTRACT 

We present narrow -band [O lll]A5007 and Ha images, as well as long-slit high-resolution 
echelle spectra of the planetary nebula NGC 6058. Our data reveal that NGC 6058 is a multi- 
polar planetary nebula of about ~ 45 arcsec in extent and formed by four bipolar outflows that 
are oriented at different position angles. Assuming homologous expansion for all the struc- 
tures, and a distance of 3.5 kpc, we obtain polar velocities around ~ 68 km s^ for three of 
them. The estimated kinematical ages suggest that the three oldest outflows have been ejected 
in intervals of ^ 1100 and ^ 400 yr during which, the ejection axis has changes its orientation 
by ~ 60° and ^ 40°, respectively. Although a inner ring-like structure is suggested by the 
direct images, its kinematics shows that no equatorial ring or toroid exists in the nebula. At the 
contrary, the long-slit spectra reveal that the ring-like structure corresponds to a fourth outflow 
that is oriented almost perpendicular to the other three. This fourth outflow is the youngest one 
and appears to be interacting with the other three, creating a protruding zone that sweeps ma- 
terial in a region almost perpendicular to the major axes of the oldest outflows. This structure 
also presents two bright arcuate regions along the same direction of the older outflows, and at 
opposite sides from the central star. From our model, we suggest that NGC 6058 could be an 
intermediate evolutionary stage between starfish planetary nebulae and multipolar planetary 
nebula with apparent equatorial lobes. 
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1 INTRODUCTION 

A Planetary Nebula (PN) is a physical system composed by a 
gaseous shell surrounding an evolved star. The shell is a prod- 
uct of the evolution of the star which is often referred as 'nu- 
cleus' or 'central star of the PN'. This system corresponds to a 
particular stage in the evolution of low- to medium-mass stars 
(M ^ 8M0), between the Asymptotic Giant Branch (AGB) 
and the white dwarf phases. The interacting stellar wind model 
(Kwok, P urton & Fitzgerald Il978r) and its generalization jBalicld 
Il987. : .lcke et alj|l989f) account for the simplest PN morphologies, 
namely, spherical, elliptical, and bipolar. However, it is well known 
that PNe show more complex morphologies, including, e.g., jets, 
multiple outflows, and point-symmetric structures (see. Acker et al. 
ri992l : ISchwarz. Corradi & Melnickll 19921: iManchado et al.lll996ah . 



In some cases, extremely complex morphologies have been found 
in images obtained with the Hubble Space Telescopqj. Thus, in the 
last decades, a renewed interest has risen to explain these complex 
morphologies. 

Among these co mplex morphologies, 

Manchado, Stanghellini & Guerrero! Il996bj defined a new 
moiphological class of PNe they called 'quadrupolar', based on 
the presence of two systems of bipolar lobes in the same PN, 
which present different orientations but share the same center. 
Examples of quadru polar PNe include K3-24, Ml-75, M2-46 , 
M3-28, and M4-14 J Manchado. Stanghellini & GueiTero"l996bj); 
NGC 6881 JGuerrero & Manchado 1 998 ; Kwok & Su 2005| 
PNG126.6+01.3 JMampaso et alj [2003) . Kn26 JGuerrero et al.l 



* E-mail:fguillen@astro.unam.mx 
© RAS 



^ Many examples of PNe and other objects can be seen at 
[http://heritage.stsci.edu/gallery/gallery.html 'and related sites. 
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Figure 1. Greyscale displays of NGC 6058. Upper (lower) panels show Ha ([O III]) images with arbitrary contrast chosen to enhance the different features 
present in the nebula. 



|2013 |). and NGC 6309 d Vazquez et ^11200 8*). More extreme is the 
case of PNe that show more than two bipolar outfl ows, the so 



called multipolar or polyp olar PNe (e. g. NGC 2440 [ Lopez et al 



199^, |V 



/azguez et • 



aiTfTsigft 



200C ; J 320 iHarman et alJ 



He 2-113, 



Sahai^ Nvman. & Wotter 



2004 NGC 5189. isabin et alj bolji 



Different explanations have been used in literature to explain the 
presence of more than one system of bipola r lobes, as, for instance , 
multiple ejections from a binary star (e.g. IVelazguez et alj|2012h . 
and bipolar ejections collirnated b y a warped irradiated disk (e.g. 
iRiikhorst. Mellema. & Ickg 120051) . In order to impose constrains 
on the models for the formation of complex PNe, the analysis of 
their spatio- kinematical structure has prove n to be cru cial (e.g. 



Balicket al. 



[1987; Miranda & SqiJ_ 



1993l;lL6pezet al. 1998: Miranda et al . 



19921: ICorradi & Schwar: 



u 



200^; 'Guerrero et al.'200: 



Vazquez et all 2008: Vav tet et alj |2009: Garcia-Diaz et al. 2009|: 



Contreras et al. 2010). In particular, high-resolution, long-slit spec- 
tra combined with narrow-band images is a powerful technique to 
obtain the 3D structure of PNe. 

In this work, we analyze the morphokinematical struc- 
ture of NGC 6058 (PNG064.6-I-48.2), located at RA(2000) = 
16''04™26':5, Dec. (2000) = +40°40'56", which is a 
relatively small PN with an apparent e lliptical morphology 
dManchado. Stanghellini & Guerrerol Il996a), but also classified 
as a multiple-shell PN by Chu et alj ( 119871) . Based on a kine- 
matical study, Sabbadiij i ll 9841) concluded that NGC 6058 is 
an evolved PN in the optically thin phase. He reports expan- 
sion velocities of 26kms~^ in [Olll] and 23kms~'^ in Ha. 



NGC 6058 is a high-excitation PN JGurzadvan & EgikvanI [19911) 
with relatively low electr on densities (Afc[On]~ 1400 cm~^, and 
iVc[ArIV] ~ 912 cm~^. I\yang et aljr2004i) and electron tempera- 
ture To ~13,200±500K JKaleilll985l) . Although NGC 6058 ap- 
pears in a relatively large number of articles (see referenc es in SIM- 
BAD), there is not a detailed analysis of its internal kinematics and 
morphology that allows us to deduce its 3D structure and the phys- 
ical processes involved on its formation. 

In this paper we present a morphokinematical study of 
NGC 6058, based on deep imaging and high-resolution, long-slit 
spectroscopy. Our study reveals that NGC 6058 belongs to the 
group of multipolar PNe mentioned above. In Section 2 we present 
our observations and results which are discussed in SectiorjS] 



2 OBSERVATIONS AND RESULTS 

2.1 CCD direct images 

Direct narrow-band images of NGC 6058 were obtained on 2007 
April 18 with the 1.5-m Harold Johnson telescope at the San 
Pedro Martir Observatory (OAN-SPlvfl). The detector was a 

■^ The Observatorio Astronomico Nacional (OAN-SPM) is located at the 
Sierra de San Pedro Maitir, Baja California, and is operated by the In- 
stitute de Astronomia of the Universidad Nacional Autonoma de Mexico 

(UNAM). 
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Figure 2. Greyscale Ha (left) and [O III] (right) images of NGC 6058 after being processed witli an uns harp-masking filter. High spatial frequencies are 
remarked with this technique. The pairs of arcs mentioned in the text are labeled as El, E2, E3, and E4. 



2048x2048 Marconi CCD (13.5 ^m pixel size) with a plate scale 
of 0. 14 arcsec pix^ ^ . The detector was set to binning 2x2, resulting 
in a plate scale of 0.28 arcsec pix^^ and a field of view of 4! 8x4.' 8. 
Ha (Ac = 6563A, FWHM = lOA) and [Olll] (Ac = 5007A, 
FWHM — 50A) filters were employed to acquire the images. See- 
ing was around 2 arcsec during observations. Total exposure times 
were 2100 s and 3600 s for the Ha and [O III] images, respectively. 
Images were processed using standard techniques of IRAF and they 
are shown in Fig.[T] 

In the Ha and [Olll] images, the nebula resembles an el- 
lipse (42 x 48 arcsec^), a nd recalls the structure previously re- 
ported by Ichu et al.l digST ) (internal shell), characterized by faint 
polar regions and a bright elliptical ring-like region. However, a de- 
tailed inspection of the images (Figure [TJ reveals the existence of 
two point- symmetric pairs of faint arcs-like filaments beyond the 
brighter emission (panels a and d), which have not been reported 
before. 

In order to enhance these structures, we have obtained 
unsharp-masking images in Ha and [O III] that are shown in Fig.|2] 
These images reveal three pairs of point- symmetric arc-like fila- 
ments, labelled El, E2, and E3, and the ring-like region, labelled 
E4. The pairs El and E3 correspond to the new point-symmetric 
arcs detected in the Ha and [O III] images, while E2 seems to trac e 
the polar caps of the elliptical shell found by Chu et al.J ( ll987D . 
Thus, El, E2, and E3 seem to be the polar caps of ellipses whose 
semi-major axes have different positions angles. The pair El is ori- 
ented at PA ~ 7° and presents an extent of ~ 38 arcsec. E2 is ori- 
ented at PA ~ —32° and its extent is ~ 33 arcsec. Between El and 
E2 two enhanced emission regions (denoted Kl and K2, Fig.[2ll can 
be distinguished at opposite sides of the central star. The pair E3 is 
oriented at PA ~ —52° and presents an extent of ~ 44 arcsec. The 



ring-like structure E4 presents a size of ~ 16 arcsec x 27 arcsec, 
with the major axis oriented at PA ~ 65°. Contrary to the other 
pairs, the polar caps of E4 are only marginally detected, whereas 
the regions around its minor axis are the brightest ones of the neb- 
ula. 

The images strongly suggest that NGC 6058 is a multipolar 
PN containing three bipolar outflows at different orientations and a 
equatorial ring-like structure that is tilted with respect to the line of 
sight. In some respects, the new images of NGC 6058 show a struc- 
ture that is reminiscent of starfish PNe (Sahai 2000), that present 
multipolar lobes accompanied by an equatorial ring-like structure. 
However, as we will see below, the kinematics of NGC 6058 shows 
that important projection effects are involved in this PN and that 
E4, rather than a equatorial ring, corresponds to a limb-brightening 
barrel-like shell, simil ar to that seen in the Ha image of NGC 7354 
dContreras et al.ll2010l) . 

We have also obtained an [O IIl]/Ha image ratio that is pre- 
sented in Fig. [3] This image shows enhancements of the [O III] 
emission at the arcs El, E2, and E3. According to Medina et al.! 
( :2009) . the [O III] enhancement indicates the presence of shocks at 
the locations of the arcs, which suggests that El, E2, and E3 are 
related to coUimated outflows. 



2.2 Long-slit echelle spectroscopy 

High-resolution, long-slit spectra were obtained o n 2003 June 6 
with the Manchester Echelle Spectrograph (MES; iMeabum et al.l 
l2003h in the 2.1-m telescope at OAN-SPM. A Site 24-^mpix~^ 
CCD with 1024 x 1024 pixels was used as detector in 2 x 2 binning 
mode (0.6 arcsec pix~^ plate scale). The spectra were centered at 
the Ha emission line using a filter ( AA=90A) to isolate the 87* or- 
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der (0.1 Apix ^ spectral scale). Exposure time was 900 s for each 
spectrum. 

A second set of spectra were obtained on 2007 July 17 and 18 
using the same telescope, instrument and detector described above. 
In this case, the spectra were centered at the [O III]A5007 emis- 
sion line using a filter ( AA=50A) to isolate the 1 W"" order (0.08 
Apix^^ spectral scale). In this case, exposure time was 1800 s for 
each spectrum. 

Data were calibrated using standard techniques for long-slit 
spectroscopy of IRAF. The slit width was 150 /xm (1.9 arcsec). The 
resulting spectral resolution (FWHM) is ~ 12kms~^ (accuracy 
±lkms~^), as measured from the lines of the ThAr calibration 
lamp. Seeing was around 1.5 arcsec during the observations. 

Figure |4] shows the slit positions used to acquire the spectra, 
superimposed on the [O III] unsharp-masking image and labeled as 
follows: SI (PA = 7°), S2 (PA = -23°), S3 (PA = -55°), and 

54 (PA = 72° ), for the spectra taken in [O III] (all of them crossing 
the central star), and S5, S6, and S7 for the spectra taken in Ha, all 
of them oriented at PA = 90° with S6 crossing the central star and 

55 and S7 shifted 7 arcsec towards the North and the South, respec- 
tively. The corresponding Position- Velocity (PV) maps are shown 
in Figures|5]and|6]for [O III] and Ha, respectively. In Figure|5] rows 
in the PV maps array correspond to the slit positions SI, S2, S3, and 
S4. Columns I to 3 show the same spectrum with different contrast 
levels. Velocity axes of all the PV maps are relative to the systemic 
velocity whereas spatial axes are relative to the location of the cen- 
tral star, or its Right Ascension, in the cases of S5 and S7. The 
systemic velocity was obtained from the velocity splitting of the 
spectral lines observed at the position of the central star. We obtain 
Vlsr — 23± 1km s^^ (VhcI — 6±lkms^^), in good agreement 
with t he value of ~ 20 ± 3.3kms~^ reported bv lSchneider et al.l 
( 119831) . In addition, the mean value of the expansion velocity, mea- 
sured at the position of the central star, is Vcx =37 km s~ ^ . 

On the PV maps (Figure |5}, the arc-like structures El, E2, 
and E3 identified in the images (Fig.|2j right) have a corresponding 
structure in the spectra. Each pair appear as the tips of an appar- 
ent velocity ellipse. In the case of El, emission can be identified at 
SI with an extent of ±19 arcsec from the central star. The radial 
velocities at the tips of El differ by 7.6knis~^, with the northern 
region being redshifted and the southern region blueshifted. E2 can 
be identified at S2 and S3 at ±16.4 arcsec from the central star. 
The northwestern tip of E2 is blueshifted by 4.2 km s^^ while the 
southeastern tip is redshifted by the same amount. We note the exis- 
tence of two bright knots on the PV map separated by ~ 33 arcsec 
and at the systemic velocity, which correspond to the knots Kl and 
K2 (see Fig.[2ll. Finally, E3 is identified at S3 as two extended fea- 
tures with a spatial extent of ±23 arcsec and a difference of radial 
velocities of 36 km s~^ . The northwestern region of E3 redshifted 
whereas its southeastern region is blueshifted. 

The PV maps at SI, S2, and S3 also show bright "knots" lo- 
cated close to the systemic velocity. These bright "knots" corre- 
sponds to the bright arcs of the ellipse identified in the image as 
E4. Remarkably, the radial velocity of these features is very similar 
to the systemic velocity. This indicates that the bright arcs should 
move mainly in the plane of the sky and, therefore, it rules out 
that E4 corresponds to a tilted ring-like structure. Moreover, the 
bright "knots" appear on the PV maps as a part of a velocity el- 
lipse, indicating that they belong to a three dimensional structure. 
The velocity ellipse is well defined at S 1 and S2 but presents strong 
distortions at S3 and S4. In particular at S4, the velocity ellipse 
seems to be open and its apparent polar regions present a complex 
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RELATIVE RA (arcsec) 

Figure 3. Greyscale image of the [O III]/Hq ratio of NGC 6058. High and 
low values of this ratio are in white and dark gray, respectively. 



kinematics. These regions correspond to the polar regions of E4 in 
the images, which are noticeably fainter than the arcs. 

Thus, the PV maps reveal the existence of a distorted probably 
ellipsoidal structure that can be only identified in the PV maps but 
not in the images. This structure, which we continue labeling E4, 
is tilted in such a way that its NE part is blueshifted and its SW 
part is redshifted. Maximum radial velocities of ~ ±53kms~^ 
are observed at S4 where a maximum angular size of 24 arcsec 
is also observed. This implies that the major axis of E4 should be 
oriented close to PA 72°, where the maximum distortions are also 
observed. This orientation is almost perpendicular to that of the 
point-symmetric arcs El, E2, E3. 

The Ha emission feature at S6 (Fig.|6j shows a structure very 
similar to that observed in [O III] at S4 (Fig.|5}. In particular, a tilted 
velocity ellipse is observed with the eastern region blueshifted and 
the western region redshifted. In addition, a marginal evidence of 
the expansion of the East and West sides of E2 is observed in S5 
and S7 slit positions, respectively (Vcxp — 30kms~^). 

A final note about the spectra around Ha is that the high- 
excitation nature of NGC 6058 results evident in them, given that 
there is not any detection of [N II]A6548, 6583 emission lines 
in these spectra. On the contrary, the high-excitation He IIA6560 
emission line is marginally detected in the spectra, with an inten- 
sity of ^ 3% that of Ha. This implies that the contribution of 
He IIA6560 to the Ha image (Fig.lTJ can be considered negligible. 



3 DISCUSSION 

From the morphological and kinematical data mentioned in the last 
section, we have identified three bipolar outflows (El, E2, and E3) 
whose small differences between the radial velocities of their polar 
caps suggest that these outflows were ejected almost in the plane 
of the sky. A fourth structure (E4), that easily can be misidenti- 
fied as an equatorial ring in the images, is clearly recognized as an 
ellipsoidal-like shell in the PV maps. The corresponding features 
of this ellipsoid in the PV maps show a higher difference in its 
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Figure 5. Position- Velocity (PV) maps of the [O III]A5007 emission line corresponding to the SI, S2, S3, and S4 slit positions (rows). Three different grey 
levels (columns) were arbitrary chosen in order to highlight different features in the maps. Kinematical features corresponding to El, E2, E3, and E4 are 
labelled in the figure. 
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Figure 4. Unsharp-masked image of NGC 6058 in [O III]. Lines indicate the 
different slit positions tliat were used to obtain our high-dispersion spectra. 
Slits labeled as SI, S2, S3, and S4, correspond to the spectra centered on 
[OlII]A5007, whereas slits labeled as S5, S6, and S7, correspond to the 
spectra centered on Ha. 



maximum radial velocities, suggesting that its orientation is very 
different with respect the outflows El, E2, and E3. 

Based on the apparent elliptical shape of the outflows in both, 
images and PV maps, we have modeled the structures El, E2, E3, 
and E4 as ellipsoids. As expected for these structures, we have as- 
sumed that each ellipsoid projects an ellipse in the plane of the sky 
which is fitted using IRAF-I-DS9 software. The ellipses were ob- 
tained trying to get the best size and shape fits for the polar caps 
El, E2, and E3, and for the bright arcs of E4. Figure |7] shows the 
fits for the case of the direct image. Using these ellipses as the pro- 
jection in 3D of the ellipsoids, we fit other parameters (inclination 
angle, velocity field, etc.) in order to reproduce the PV maps. 

For each ellipsoid, the apparent polar radius and velocity are 
measured from the image and the corresponding long-slit spectrum, 
whereas the equatorial radius and velocity are obtained from the 
fitted ellipse and the kinematical model, respectively. For the model 
we have always assumed a homologous expansion velocity {v oc r) 
and small deviations (< 20°) around 90° for the inclination angle 
(i). A sketch showing the main parameters of an ellipsoidal shell 
(a), as well as the observer view for the image (b) and the PV map 
(c), are presented in Figure |8] With the considerations mentioned 
above, and following the nomenclature in Figure Ula), it can be 
demonstrated that 



tan a — 






where r^ is the equatorial radius (semi-minor axis); Vp^ is the ap- 
parent polar radius (projected semi-major axis); v^ is the equatorial 
velocity, and Vpx is the apparent polar velocity (radial component). 
The inclination angle of the ellipsoid major axis, with respect to the 
line-of-sight, i, can be derived as i = 90° + a, if the North pole of 
the ellipsoid is moving away from the observer, and as i = 90° — a 
if coming closer. 

It is important to mention that, as can be noted in Fig [8] the 
radial velocity measured at the position of the central star, Vcx , is 
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Figure 6. Position- Velocity maps in Ha. Labels correspond to slit position 
shown in Fig.|4] East is up and West is down. 



different than Ve as well as Vex- However, using a rough approxi- 
mation of the polar-to-equatorial ratio r-p/ro ~ 2 (which is reason- 
able, given the semi-axes ratios of the ellipses), assuming a "small" 
value for a ~ 20°, and using the cartesian equation for a vertical 
ellipse, we have estimated the ratio 



Wo 



rp/re 



1.047 



''c ysin^ a + {rp/vc)^ cos^ a 

i.e., an error < 5% is commited when considering the approxima- 
tion Va ~ Vex- Therefore, we will use this approximation and its 
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Figure 7. We present the four ellipses model that best fit the morphology 
observed in NGC, 6058. The ellipses are superimposed on the [O III]A5007 
unsharp-masking image shown in Fig.|2}). 



associated error to estimate a, i, and the other parameters in the 
following calculations. 

Polar radius and velocity can be now obtained using the re- 
lationships for the projected Vp and the homologous velocity law, 
respectively 



:; Vp = Vc { — ] = v^x \ 2^ j , 

I V Tc / V re sm 1 1 



In addition, the kinematical age r^ can be derived as 



4744 X 



r '^p 1 




r D 1 

.kpc. 




r "p 1 


. arcsec . 


.kms~^ . 



where D is the distance from the Earth to the nebula (in kpc). Val- 



ues of i, 



, and Tk , obtained for all the ellipsoids and assuming 



a distance of 3.5 kpc dCahn. Kaler & Stanghellinill 19920 . are pre- 
sented in Table [T] It is important to mention that E4 was excluded 
from Table [T] given that is not possible to obtain a suitable kine- 
matic model for this ellipse. As we will explain below, we interpret 
this as the result of outflows interaction. However, as a reference 
only, we propose the following limits for this ratio 

1.3 < IE <2. 

re 

where the lower limit corresponds to the apparent axes ratio 
(rpz/ro) seen in E4, and the upper limit was choose as a reason- 
able value given the apparent polar-to-equatorial ratios from the 
other ellipses. Using these values, the possible inclination angle for 
such ellipse would be between 50 and 130° . Consequently, a kine- 
matical age for E4 can be estimated as Tk ~ 3400 yr. However, as 
we mentioned above, this value must be taken cautiously. 

As can be noted in Table[T] given the uncertainties in the kine- 
matical ages, it is not simple to describe the sequence of events 
that have occurred in the formation of NGC 6058. Nevertheless, 
based on the kinematical ages as well as on the orientation of the 
structures, we propose a scenario in which a series of bipolar out- 
flows could shape the morphology of NGC 6058. In our proposal, 
the first outflow was E3 along PA ~ —52°. The second ejection 



was El along PA ~ 7°, involving a change in the orientation of the 
ejection axis of ~ 60° with respect to that of El . The third ejection 
was E2 along PA ~ —30° that involves a change of orientation of 
~ 40° with respect to El. We note that E3 and El are ejected with 
their northern/northwestern half pointing away from the observer 
whereas the northwestern half of E2 points towards the observer. 

Moreover, the deduced kinematical ages (Table [T) suggest a 
time span of ~ 1100 yr between the two first coUimated ejections 
and only ~ 400 yr between the second and the third. These results 
point to the existence of episodic ejections with changes in the ori- 
entation of the main axis, suggesting precession of the collimation 
axis. About > 1000 yr after ejection of E2, a new shell is formed 
which is identified with E4. In this case, the main ejection axis has 
largely changed with respect to that of EI-E3, as the polar axis 
of E4 is almost perpendicular to the previous outflows. Our data 
strongly suggest that E4 interacts with the previous outflows. The 
bright arc-like regions observed in the direct images are probably 
a result of the interaction of the equatorial regions of E4 with parts 
of E1-E3. Furthermore, the large kinematical distortions observed 
in E4 along PA around 72° suggest that E4 also interacts with the 
equatorial regions of the previous outflows and protrudes. We note 
that this interaction may result in a deceleration of E4 and, there- 
fore, the deduced kinematical age of 3400 yr is probably an upper 
limit, implying that E4 has been formed noticeably later than the 
previous outflows. 

As a result of this study we can conclude that NGC 6058 be- 
longs to the group of multipolar PNe. It is interesting to compare 
NGC 6058 with other multipolar PNe. In particular, we find re- 
marka ble si milarities between NGC 6058 and the so called 'starfish 
PNe' JSaha i 2000). These PNe present three or four pairs of bipo- 
lar lobes and an bright elliptical ring-like structure that can be in- 
terpreted as a tilted equatorial ring. We may relate the El, E2, and 
E3 outflows with the multiple pairs of lobes of starfish PNe, and 
the bright arcs corresponding to E4 with the equatorial ring of the 
starfish. However, in NGC 6058 our kinematical data reveal that 
the bright arcs do not trace a tilted equatorial ring but a part of an 
ellipsoid shell. In this respect, it would be interesting to study the 
internal kinematical of starfish PNe in order to check whether the 
observed elliptical rings have the kinematics expected from a tilted 
ring or, at the contrary, they trace an ellipsoidal shell, as in the case 
of NGC 6058. 

It is interesting to speculate about the possible future evo- 
lution of NGC 6058. In particular, the polar regions of E4 could 
break and expand at a higher velocity, creating, after some time, a 
pair of large bipolar structures. If so, NGC 6058 would consist in a 
pair of extended bipolar structures, accompanied of smaller bipo- 
lar lobes around its equatorial region, which would correspond to 
EI-E3. This scenario could provide an explanation for some PNe 
like Hb5 and K3 -17 JMiranda. Ramos-Larios & Guerrerol l2O10t 
iLopez et al.ll2012h in which large bipolar lobes are observed ac- 
companied by secondary smaller lobes around the equator of the 
largest ones. By all these comments, we suggest that NGC 6058 
could be an intermediate evolutionary stage between starfish PNe 
and PNe with equatorial bipolar lobes. 



4 CONCLUSIONS 

We have carried out a morphokinematic analysis of the planetary 
nebula NGC 6058, based on direct imaging and high-dispersion 
spectroscopy. NGC 6058 can be described as a multipolar PN, 
formed by three bipolar outflows with different orientations and 
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Table 1. Parameters of ellipsoids measured and derived from models: position angle (PA), semi-minor axis (re), apparent semi-major axis (rp), equatorial and 
polar apparent velocities (Vc) and (Dp), inclination angle (i), semi-major axis (rp), polar velocity {Vp), and kinematical age (tj^). 





PA 

n 


(arcsec) 


rpz 
(arcsec) 


Vex 

(kms-1) 


Vpx 

(kms-l) 


i 

n 


^P 
(arcsec) 


Dp 
(kms-i) 


(yr) 


El 
E2 
E3 


+1 

-30 

-52 


11 ±1 
10 ±1 
13 ±1 


19 ±1 
18 ±1 
23 ±1 


37 ±2 
37 ±2 
37 ±2 


3.8 ±1 

4.2 ±1 
18 ±1 


93 ±2 
86 ±2 
106 ±2 


19 ±1 

18 ±1 
24 ±1 


67 ±8 
69 ±9 
67 ±10 


4 800 ± 600 

4 400 ±600 

5 900 ± 600 






Figure 8. (a) A sketch of an expanding ellipsoidal structure, whose semi- 
major axis has an inclination angle a respect to the plane of the sky, seen 
from an observer on Earth. Main physical parameters are labeled, (b) The 
elUpse (image) seen by the observer (c) The position-velocity map (spec- 
trum) from a slit crossing the major axis. 



polar velocities around ~ 68kms~^, as well as a distorted ellip- 
soidal inner shell with its major axis almost perpendicular to those 
of the other three outflows. Assuming homologous expansion for 
all the structures and a distance of 3.5 kpc, we obtain kinematical 
ages ranged from ~ 3400 to ~ 5900 yr for the different outflows. 
Our data suggest a certain systematic behaviour in the changes of 
the orientation of the main axis, which points out to episodic, pre- 
cessing ejections in NGC6058. The ellipsoidal inner shell is re- 
leased much after than the three collimated outflows and its kine- 
matics reveals clear signs of interaction with the previous ejections. 
A comparison with other multipolar PNe suggests that NGC6058 
could be an intermediate stage between (young) starfish PNe and 
more evolved PNe that show a large pair of bipolar lobes and much 
smaller equatorial lobes. 
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